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. We consider the rate of fermion-antifermion pair creation by an external field. We derive a 
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rate formula that is valid for a coupling with arbitrary vector and axial vector components to 
nuclear matter, a problem with astrophysical relevance. We present an estimate for the creation 
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rate per unit volume, compare this to previous results and comment on the role of the neutrino mass. 
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^ . I. INTRODUCTION 

a: 

(-} • Starting with Schwinger's classical account 

electric field, fermion pair creation has been the subject of continued interest. A variety of pair 



of electron-positron pair creation by an external 



creation rates for specific external fields in quantum electrodynamics can be found in the literature, 

true quantum field theory phenomenon: the creation of particles from the vacuum. 

Because neutrinos carry weak charge, one expects that an external Z-boson field can produce 
neutrino-antineutrino pairs in a similar manner. The concept of an external Z-boson field can be 



seen as arising from a distribution of nuclear matter (in the sense of ref. 



ll|). Neutron stars are 



a prime example of such a matter distribution and their neutrino emission by this mechanism was 
studied using non-perturbative methods [ill I12I ll.j . Pair creation of neutrinos is also studied 
in relation to the stability of neutron stars, see Ref. ^J] and references therein. Although Refs. 
[ill I12I I13I ] find typical neutrino fluxes that are too small to be observable, we believe it is 
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worthwhile to study such a relatively unexplored neutrino source from a different point of view. 
In particular, we want to develop a method that is not limited to a specific source but allows us 
to draw conclusions with a broad applicability. This can then be applied to study e.g. neutrino 
pair creation by non-standard model weakly interacting particles or domain walls. 



In the present letter, we study the creation of neutrino pairs in a perturbative way. We present a 
first order computation of the pair creation rate per volume, with a dynamical nuclear configuration 
acting as a source. The reasons for using perturbation theory are twofold. First, the axial coupling 
to the Z-boson prevents an easy generalization of non-perturbative QED methods. Second, non- 
perturbative methods generally consider a very specific source, or class of sources, from the start. 
The perturbative method is more flexible in the sense that a specific source is folded in at the 
end. This allows us to keep separate the physics of the pair creation process and that of a specific 
source. 

In part, our computation was triggered by the results presented in ref. 11], in which the 
creation of neutrinos by a time-dependent nuclear distribution is studied. One of the results 



in ref. 



is that the overall rate is proportional to the square of the neutrino mass. This 
implies that there can be no pair creation of massless neutrinos. The question arises whether 
this is a manifestation of a general chiral suppression mechanism or a consequence of the specific 
source considered. We shall see that the perturbative viewpoint contributes to a more complete 
understanding of this effect. 

The paper is organized as follows. In section [H] we discuss the theoretical background of pair 
creation processes for fermions and introduce the relevant quantities. In section ITTll we discuss the 
perturbative computation. The result is then applied to neutrinos in section IIVI and we present 
our conclusions in section Ivl 



II. PAIR CREATION PHYSICS 



We study fermions that are coupled to an external source j. The interaction Lagrangian reads 

C irA =j„{x)^{x)T^{x). (1) 



The source is fully prescribed and has no further dynamics. We choose the coupling of the general 
form 

r" = r(c v -c A1 5 ), (2) 

where cy (ca) is the vector (axial vector) coefficient; the coupling constant is absorbed in j. 



Following ref. [Ijj, we introduce the overlap of asymptotic 'in' and 'out' vacua to describe the pair 
creation process: 

S (j) = (0, oo|0, -co), - (0, oo|S|0, oo)j , (3) 

where S is the scattering operator and the subscript is a reminder that a source is switched on 
and off adiabatically somewhere between t = — oo and t = oo. The probability that a system that 
started in the vacuum state will remain in the vacuum state is then expressed |i!t| as: 

|(0,oo|0,-oo) J -| 2 = exp(-W r )=exp(- (d^xwix)). (4) 



For a positive W , this probability is between zero and one which signals a non-zero probability for 
the creation of a fermion pair. Now suppose that w(x) = w is constant. We can embed the system 
in a box of size V x T, write W = wVT and choose the box small enough such that W < 1: 

|(0,oo|0,-oo)j| 2 ~ 1-wVT, (5) 

which supports the interpretation of the function w(x) as the probability per unit time and volume 
to create a pair at space-time location x. Such a rate density is the physical quantity of interest. 
For QED, the Schwinger formula [| states that for a photon field of the form A^(x) = j^(x) 
Ml 0.1). (I 'I;. 
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where m is the electron mass. We mention that refs. 



HQ 



conclude that this result extends to 



the case of neutrino pair creation by a source of the same form. 



To compute the rate density, we use perturbative quantum field theory: 

(0 1 oo|0,-oo) J - = Z[j] = eotp (iW\j]), 



(7) 
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where W[j] is the generating functional of connected n-point functions. 1 In this context, W[j] is 
also the effective action for the external field j [lfl . 

The interaction Lagrangian Q only contains a vertex that couples to the external field. 
Therefore W[j] represents an infinite sum of fermion loop diagrams, labeled by the number of 
vertices which are all connected to the external field. In terms of W that was introduced in eq. 

m, 

W = 2lmW\j}. (8) 

The fermion loop diagram with one external field vertex is zero by momentum conservation, so 
the first non-zero contribution is from the loop with two external field vertices, i.e. the two-point 
function. This is the object that we will compute in section lllTl Its contribution to the pair creation 
rate W2 is found by folding in the sources according to the formula 2 

W 2 = - J J^L^^-p) TmJT'ip) , (9) 

where represents the two-point function, with prefactors as chosen in eq. Hl()|) . For time-like 
currents, j l i{p)jv{—p) ImS M "(p) < since a probability cannot exceed one. For a given W, the 
pair creation density follows by extracting the function w(x). 



There has to be enough energy in the source to put two virtual particles on-shell. For the pertur- 
bative mechanism that we describe, this implies a threshold energy for the source insertions. This 
is in contrast to the non-perturbative effect, which can be thought of as an infinite sum of loop 
diagrams with an increasing number of source insertions. This infinite amount of sources conspire 
to create a pair and the amount of energy per source insertion can be arbitrarily small. 

For QED it is known that the real part of the sum of loop diagrams has a divergent structure, 
which can be used to extract non-perturbative results by performing a Borel transformation |lfj . 
We do not know whether or not a similar procedure can be applied in this more general situation. 



1 The use of W and W[j] may be confusing, but both symbols are standard in the literature. The generating 
functional will always be denoted with its argument j. 

2 We use a metric tensor g^ v = diag (1, —1, —1, —1) throughout this paper. 



k+p 



FIG. 1: Fermion loop diagram with two external sources attached. The external field couples directly 
(i.e. without propagators) to the loop. 



III. THE TWO-POINT FUNCTION 

The two-point function without external sources is transcribed from figure ^ We find that, in 
dimensional regularization with n = 4 — e, 

n r d n k tr [(# + m) T» (jk + p + m) V] 
J (2ir) n (k z — m z + ie)((k + p) — ra l + ie) 

where m is the fermion mass. From eq. we are interested in the imaginary part of this 

expression, which is finite. Note that we integrate over the fermion momentum; in the source's 
rest frame (where the particles are created back to back), the fermion and the antifermion each 
carry half of the energy. 

Expression IjlUI) is reduced to a linear combination of scalar integrals in the fashion of 
Passarino-Veltman [17|. A series expansion in e reveals the divergent structure, and after some 
algebra the problem depends only on the one- and two-point scalar integrals. The one-point scalar 
integral is real, the two-point integral develops an imaginary part if p 2 > 4m 2 which means there 
should be enough energy in the source to create two fermions. If this is not satisfied, E^" is purely 
real and there is no pair creation. The final result is the following expression: 

IniE"» = -^- 2 [(4 - c\) Ef (p) + (4 + c\) E£(p)] Q(j? - 4m 2 ) , (11a) 



Ef&o = w t /i-^r, (nb) 



p 

2 



E£(P) = ^ (pV - P> V - - ) Jl ~ ^ • (He) 

For some typical values of the parameters cy and ca, this result can be compared to the literature 

3D 
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From expressions (|llfl we observe that for massless fermions Yf^ — 0, so that only the second term 
contributes. This means that the physics is qualitatively insensitive to different choices of cy and 
ca] only the square sum is quantitatively important. We conclude that the difference between the 
two-point functions with two different normalized sets of couplings (e.g. purely vector, purely axial 
vector) is proportional to m 2 . 

The contribution due to the three-point diagram should be interpreted with care. In QED 
it vanishes by Furry's theorem, but for axial couplings it contributes to the axial anomaly. This 
means one should verify that the final result does not depend on the regularization procedure. For 
the present calculation, this is not an issue. 



IV. NEUTRINO PAIR CREATION TO FIRST ORDER 
A. The general case 

We specialize to neutrino pair creation by putting cy — ca — 1/2 in the expression for the two-point 
function Qllfl. Combining eqs. © and i|ll|) . we find 



^ = ~ J^J(P 2 -*rn^l-^[MP,j)+rn 2 K(p,j)] , (12a) 

MP,J) = P 2 W) ■ 3(-P)} - [P ■ J(P)) (P ■ j(-p)} > (12b) 

2m 2 

m Fi(p,j) = -m [i(p)-i(-p)]--rb-iW)(P'i(-p)] • ( 12c ) 

Without loss of generality, we consider a source with a density component and a spatial current in 
the z direction: 

3n(p) = (jo(p))0,0,i 3 (p)), = (J5,pr,p 3 ). (13) 

Here p M labels the energy and momentum of the source. Though the current is directed in the z 
direction, we allow for a dependence on the transverse direction by leaving pV unspecified. The 
two terms in (I12|l can be written as 

MP,J) = -Pt (bol 2 - |j 3 | 2 ) - \Ej 3 -p 3jQ \ 2 , (14a) 
m 2 ^(p,j) = F2 _ g _ 2 (£; 2 |jo| 2 +P3 2 |J 3 | 2 ~ E P3 (j f 3 +j* j 3 )) 

Hi Pj, p 3 

-m 2 (|j | 2 -|j 3 | 2 ) • (14b) 

We do not simplify these equations further, as we do not want to constrain the properties of the 
source. 



It is instructive to analyze the massless limit in more detail. In this case only J-q in 112|) 
contributes, so that 

1 f d A p 



W 2 (m = 0) = -— j [p 2 (j(p) ■ j(- P )) - (p ■ j(p)) (p ■ j(-p))] . (15) 

In analogy with QED, we introduce a field strength F^ v (p) = ipujp,(p) — ipp,ju(p) and its 'electric' 
and 'magnetic' components Ei and Bi and find: 

W 2 (m = 0) = ~ [^(p)F^(-p)] (16a) 

= ^ /(04 (P)^("P ) - ^W^(-P)] • (16b) 



This is exactly half of the QED result |lf| if we insert a factor e 2 from the coupling constants, 
which reflects the discussion in the previous section. In electrodynamics, E and B are the physical 
electric and magnetic fields and one can go to a frame in which B — 0. Then eq. i|16|) yields a 
positive result from which we conclude that the creation of massless particles by the two-point 
mechanism is in general possible. Eq. (|16fl is consistent with the massless limit of the first-order 
effective action in an axial background that was computed in ref. jlfll |. 

It is interesting to compare this result to the creation of neutrinos by an external electromag- 
netic field as computed in ref. 20]. In that case, the pair creation rate is proportional to m 2 and 
depends on the electromagnetic invariant E ■ B. 

B. The time-dependent density 

We consider a time-dependent distribution of nuclear matter, described by the following source 
term: 3 

U*) = G F /V2(n^(l-j 5 )n) = (j (t), 0,0,0) , (17a) 

j (t) = -j=n N (t), (17b) 

where njy is the number density of the nuclear matter distribution and Gf is Fermi's constant. 
This is the specific background that we refer to as a time-dependent density. Our main motivation 
for this source is to compare the perturbative results with the non-perturbative results of ref. jll| . 



3 This source originates from an effective four-fermion description, see ref. Illl . Note that contains the axial 
current; since the neutrons are massive, axial symmetry is broken and the current need not be divergence-free. 
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For simplicity (and because any source can be decomposed into a trigonometric sum) we 
assume a monochromatic source: jo(t) = Eq cos tot. In Fourier space, this is 

Jo(p) = ^(27r) 4 S(p)[6(E-iu)+S(E + u;)] 1 (18a) 
E = %n w (0). (18b) 

Inserting the source l|18|) into eq. results in products of delta functions. We employ a box 
normalization procedure to reduce these to a single delta function and a factor V x T and find 

W 2 = -VT Gp2 ^ w)2 [lm£°V;p = 0) + Im£ 00 (-co;p = 0)1 . (19) 
8 

Using eq. (fTTJ) . with cy = ca = 1/2, we see 



ImE m (i,; P - 0) = Jl - ^ , (20) 

87T V W 

leading to the following pair creation probability per unit time and volume: 

W? m 2 I 4m 2 



^ = ^=3^V 1 -^ W - (21) 

The rate density scales with the square of the nuclear density, as expected for the two-point 
mechanism. 

We see that the rate is proportional to to 2 , which could have been anticipated from eq. 
I|14(l because the time-dependent density (jl7|l is characterized by j'3 = pr = P3 = so that 
•^o(P)i) = 0. Rcf. 11] also finds the m 2 proportionality for sources with a time-dependent current 
in the z-direction. Equation l(Ti|) suggests that such sources can contribute to first order for a 
zero neutrino mass. 

To derive an order-of-magnitude estimate for the number of created neutrinos per unit volume per 
unit time, we take the square root factor in eq. (|21() of order unity, use a neutrino mass of 0.1 eV 
and assume a 'reduced density' Gftt-n I 'V% ~ 1 eV, such as in a neutron star |ll| : 

= (0.1 eV) 2 2 _ i -4 eV 4 ~ 10 26 s" 1 cm" 3 . (22) 

At the pair creation threshold, this corresponds to an energy output of order 10 13 erg cm -3 s . 
Ref. |ll| estimates the energy output of neutrinos that are created non-perturbatively by an 
oscillating neutron star to be of order 10 3 erg cm' 3 s -1 . However, these numbers should not be 
compared because the (realistic) driving frequency that is considered in \l\ is so low that the 
perturbative mechanism is not operational. 



9 

As follows from cq. there can only be pair creation by the two-point mechanism if 

uj 2 > 4m 2 . With a neutrino mass of 0.1 eV, the creation of a neutrino-antineutrino pair requires 
a driving frequency of at least 3 • 10 14 Hz. The coherence length of such a system is roughly 10 -4 
cm, so it is not very feasible to look for an oscillating astrophysical object that would produce 
an appreciable amount of neutrinos with this mechanism. However, the value of our computation 
lies in its general applicability. We are not limited to this particular type of sources, and we 
believe it may be interesting to study sources of a more transient nature such as a forming neutron 
star. Alternatively, one could consider weakly interacting particles beyond the standard model or 
domain walls as a source. 



V. CONCLUSION 



We have described pair creation of fermions by an external field to first order in perturbation 
theory and found the contribution by the two-point mechanism for a general coupling. Our main 
result is eq. which should be interpreted in the context of eq. ©■ We observe that at this 

order in perturbation theory, the difference in pair creation rates between two sets of normalized 
coupling coefficients {cy, ca} is proportional to the square of the fermion mass. 

For the case of neutrino pair creation by a distribution of nuclear matter, we have derived 
expressions (|12|) and (|14|l . From this result we observe that, to first order, neutrino pair creation is 
possible with a suitable source if neutrinos would have been massless particles. We then considered 
pair creation of neutrinos by the time-dependent density of eq. i|17|) . For this specific source, we 
conclude that the production rate due to the two-point contribution l|21|) is proportional to the 
square of the neutrino mass. This conclusion is in qualitative agreement with the non-perturbative 
result derived in jll |. 

The method that we presented in this paper is suitable to study different types of sources 
and we intend to do so in the future. 
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